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FOREWORD 


This  is  Part  1 of  the  final  report  under  Contract  DEEA-76-C-0019 
entitled  Atmospheric  Transmittance  Study  with  the  Meteorological 
Satellite  Technical  Area  at  White  Sands  Missile  Range.  Part  II  contains 
the  report  on  the  study  of  the  inversion  of  the  radiative  transfer 
equation  for  the.  temperature  profile  as  well  as  for  the  absorber  con- 
centration in  the  15y  C'0?  band.  Also  included  there  is  the  method  used 
in  this  study  for  the  calculation  of  atmospheric  transmittance  using 
line  spectral  parameters.  Part  III  deals  with  the  study  of  the  effects 
of  clouds  on  the  inversion  techniques.  Essentially,  Part  1 deals  with 
the  study  and  development  of  band  models  for  use  in  connection  with 
techniques  for  the  calculation  of  atmospheric  transmittance  along 
slant-paths . 
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t.  INTRODUCTION 


The  Inversion  ot  tho  radiative  transtei  equ.it  ton  for  the  tempera- 
ture us  In#  satellite  radiance  measurements  requires  the  use  ot  a 
transmittance  function  for  the  atmospheric  region.  The  transmittance 
function  itself,  however,  it  it  is  to  represent  accurately  the  ef- 
fects of  the  gaseous  absorbers  it  must  he  an  explicit  function  ot  the 
temperature  sought  for  in  the  inversion.  This  implies  that  the  nu- 
merical procedures  adopted  must  include  the  capability  for  iterating 
in  the  computation  of  the  transmittance  as  close  guesses  are  obtained 
for  the  temperature.  This  would  be  a relatively  simple  modification 
if  it  were  not  for  the  i naccoss i hi l i t v of  appropriate  t ransmi t t ance 
funct ions . 

The  principal  problems  associated  with  the  transmittance  function 
are:  long  cor.tput.it  Iona  1 times,  deviations  trow  measured  data,  complex 

analytical  represent  at  ions , inaccurate  comput at ieua 1 results,  and 
limitations  of  the  conveislon  from  inhomogeneous  to  homogeneous  paths. 
In  this  report  some  ot  these  factors  are  discussed  in  connect  ion  with 
the  most  recent  methods  and  models  available  in  the  literature.  based 
on  the  work  performed  with  line-bv-ltne  data,  available  from  another 
part  of  this  of  tort.  It  is  recommended  that  a transmittance  funct ion 
in  the  form  of  a polynomial  In  the  pressure,  temperature  and  absorber 
amount  be  used  as  a homogeneous  path  model.  The  model,  however,  would 
be  developed  with  the  vertical  transmittance  tor  each  laver  in  a hot  i- 
r.ontallv  stratified  atmosphere.  The  transmittance  to  each  level  is 
then  obtained  through  mul t ipl icat ion  ot  the  transmittance  through  pre- 
vious lavors. 
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2.  TRANSMITTANCE  FOR  HOMOGENEOUS  MEDIA 


Infrared  radiation  passing  through  a nonscattering  medium  is  ab- 
sorbed bv  the  molecules  along  the  path  in  the  process  of  rotational- 
vibrational  energy  transitions.  The  law  governing  the  absorption  mech- 
anism is  that  due  to  Beer,  and  its  application  is  dependent  on  line- 
broadening mechanisms  of  the  Doppler  or  I.orentz  type.  In  spite  of  the 
fact  that  the  physical  processes,  as  wekl  as  their  mathematical  repre- 
sentations are  well  known,  it  is  still  a challenge  to  arrive  a a 
practical  and  accurate  function  for  the  transmittance  through  the 
atmospheric  region.  In  this  section  a presentation  is  made  of  two  band 
models  of  the  polynomial  type,  which  have  been  found  to  be  of  signifi- 
cance value  in  atmospheric  work. 


2.1  The  transmittance  Function 

Spectral  transmittance  through  an  absorbing  gas  is  given  bv  the 

monochromatic  form  of  Beer's  Law,  that  is 

ft' 


K ((P,T)  du ' 


(2.1) 


Tv(P,T,u)  = e 

where  v is  the  frequency,  1’  is  the  pressure,  T is  the  temperature,  u is 
the  absorber  amount  and  K is  the  absorption  coefficient.  For  homoge- 
neous media  P and  T do  not  change  along  the  path,  such  that  (2.1) 
becomes 


-Kv(P,T)u 


Tv(P,T,u)  = e 


(2.2) 


Over  a spectral  band  - Vp  - Vj  use  is  commonlv  made  of  the  average 
value  of  (2.2),  which  bv  the  first  mean-value  theorem  becomes 
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Since  T.  is  the  quantity  dealt  with  in  the  work  reported  here,  tut 
thor  reference  will  use  the  short  notation 

T\vU\T,u>  T(P,T.u)  (.\4> 

The  absorption  coefficient  K^,  it'  (2.2)  depends  on  the  nature  ot 
the  line-broadening  mechanism,  on  the  number  and  arrangement  ot  the 
lint's  and  on  the  tvpo  ot  molecule.  Spec i f lent  ion  ot  the  broadening 
as  being  due  to  either  Doppler  or  lorent.*  has  allowed  tor  the  evalua- 
tion ot  (J.'),  with  resulting  simple  analytical  expressions.  However, 
the  simplicity  of  the  expressions  is  obtained  at  the  cost  of  ovov- 
stmplication  of  the  physical  processes  involved.  These  results,  or 
"band  models",  do  not  as  a general  title  represent  the  actual  trans- 
mittance to  a degree  o!  accuracy  high  enough  to  justify  theit  unbiassed 
use.  However,  use  mav  be  made  of  their  form  in  the  obtainment  ot  more 
empirically  satisfactory  functions  of  the  variables  P,  T and  l’. 


2.2  Polynomial  in  Weak- and  Strong-line  Functions 

For  the  special  case  of  statistically  distributed  Lorent* tan 
lines  over  A^  with  a Poisson  intensity  distribution  function,  Mav or 
[11  and  Coodv  [2)  evaluated  f2.3)  obtaining 
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where  the  variables  b and  t*’  are 
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and  a Is  the  moan  half  width  normal izod  to  ci^,  d is  the  mean  line 

spacing,  S is  the  mean  lino  intensity  normalized  to  S , and  a , S 

o O O 

are  their  values  at  standard  temperature  T and  pressure  P . In  t lie 

o o 

weak-line  (i.e.,  ij>  < < 1)  and  strong-line  (i.o.,  <!•  N > 1)  limiting 
conditions,  (2.5)  results  in  the  approximations 

-tv  - exp  (2.6) 

in  which  and  I are  called,  respectively,  the  weak-and  strong-line 
models.  The  pressure  and  temperature  dependence  is  introduced  in  the 
last  three  equations  through  and  S. 

A re-arrangement  of  (2.5)  leads  to  the  geometrical  form 

(r' 

which  states  that  the  inverse  of  the  natural  logarithm  of  the  two 

approximations  are  in  quadrature  with  the  complete  transmittance 

function.  If  this  result  is  considered  universal  for  transmittance, 

then  (2.81  may  he  generalized  with  the  choice  of  a 1 which  is  valid 

also  for  regularly  distributed  lines.  The  choice  of  T as  in  (2.6) 

w 

is  a valid  selection  for  any  type  of  line  arrangements. 

Of  the  strong-line  models  available  the  most  adaptable  to  (2. 8' 


is  the  one  due  to  Kino  [3),  which  represents  a continuous  distribution 
of  Lorentzlan  lines  from  regular  to  random  through  the  variation  of  a 
parameter.  Accordingly, 


(2.9) 


n T (n)  ( 


l/n 


where  G is  the  Incomplete  gamma  function,  F is  the  regular  gumma  func- 
tion, and  n Is  the  parameter  that  specifies  the  type  of  line  distribu- 
tion (e.g.  n “ 0.5  for  regular  and  n * 1.0  for  random).  The  generali- 
zation of  the  Mayer-Goody  model  in  this  manner  is  due  to  '/.achor  [ 4 , 5 } . 

In  spite  of  the  abovi -mentioned  broad  interpretations  of  such 
classical  band  models,  their  adaptation  to  measurements  is  always 
restricted  by  hidden  factors  affecting  the  measurements.  The  deviation 
from  the  idealism  of  the  mathematical  formulation  of  the  physical 
process  of  absorption  may  be  accounted  for  by  the  addition  of  terms  in 
(2.8).  Thus,  a polynomial  model  may  be  proposed  in  the  form  [b]. 


-1 

In  x (P.T.u) 


Cj  x^  + C2  y + C3  xy  + C^x  y + C^xy 


+ ...  (2.10) 
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-l 
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TS 

and  the  C . 's  (i  * 1,.  . . , N)  are  determined  from  least-squares  curve- 
fitting to  measured  data.  The  N-term  polynomial  in  (2.10)  is  described 
by  N + 3 spectral  parameters,  two  of  which  arise  from  in  (2.9)  and 
one  from  r in  (2.6). 

The  dependence  of  t on  P,T  and  u may  he  explicitly  made  apparent 
with  the  use  of  the  normalized  Lorentzian  half-width  a in  6 and 


that  is 


6 


h 


(2. 11) 


Also,  the  model  variables  in  (2.6)  and  (2.9)  mav  be  separated  from 
their  spectral  dependence  by  expanding  in  the  form 


fty  = k Su  (2.12) 

= c sau  (2.13) 


where  k = So/d  and  C = 2 d^S^d  are  spectral  parameters.  The  complete 
set  of  spectral  parameters  describing  an  N-term  expansion  of  the  poly- 
nomial model  in  (2.10)  consists  of  C. , C„,  C,,...,  C„,  k,  n and  C.  In 

I z j N 

the  development  of  the  model  with  transmittance  data,  k,  n and  C are 
determined  first  and,  then,  the  's  are  determined  simultaneously. 

Details  of  the  process  will  be  given  in  Section  4 of  this  report. 

Although  the  model  has  been  derived  for  homogeneous  paths,  it  may  he 
used  for  inhomogeneous  paths  with  the  use  of  variables  equivalent  to 
Su  and  Sau.  These  may  be  called  the  weak-line  and  strong-line  vari- 
ables, respectively.  Details  on  these  equivalences  will  be  discussed 
in  Section  6.2. 

| 

2.3  Polynomial  in  Pressure,  Temperature  & Concentration 

j 

A totally  different  approach  may  be  followed  in  the  obtainment 
of  a polynomial  for  T in  terms  of  functions  explicitly  showing  the 
dependence  on  the  pressure,  temperature  and  gas  amount  (or  concentra- 
tion). With  the  use  of  (2.12)  and  (2.13),  the  weak-and  strong-line 
models  in  (2.6)  and  (2.7),  respectively  become 

I ■=  exp  {-k  Su)  (2.14) 
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( 2 . IS) 


r J . .ct/2  1/2  1/2  ( 

« *xp  V k S a U > (2.  IS) 

where  k‘  - *C/2.  An  Inspection  of  the  powers  associated  with  the 
variables  in  these  two  limiting  equations  for  transmi ttnncc  suggests 
the  assumption  of  a general  transmittance  equation  of  the  form 

T (P,  T,  u)  « exp  ^-k"  S'*  a*’  u‘^  (2.16) 

where  k"  , a,  b,  ami  c are  spectral  constants.  Replacing  a with  (2.11), 
using  the  approximation  (7] 

kd 


s ** 


(2.17) 


and  taking  the  natural  logarithm  of  (2.16)  twice  leads  to 

<’n  |^tn  t (P,  T,  A - Sj  + a , fn  u + a.}  (n  (P/lM  + a^  fn  (T/T^)  (2.18) 

where  d,  a ^ , a^,,  a.^,  and  a^  are  constants  related  simply  to  the 
earlier  k",  a,  h and  c. 

The  polynomial  in  (2.18)  may  be  further  expanded  to  allow  for 
better  flexibility  in  curve-f itting  to  experimental  data.  With  the 
def init ions 

x ■ fn  u 

v - Cn  (p/p  ) 
o 

z - (n  (T/T  ) 
o 

(2.18)  may  be  expanded  into  the  form 

<n  <-fn  r (P,T,u)»  • a.  + ax  + a v + a z + a xv  4 a vz  + a,zx  4 ...  (2.19) 

\ I I 2 3 s S 6 7 

which  is  a form  of  the  model  proposed  bv  Smith  ( 8 1 . 
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3.  TRANSMITTANCK  FOR  lNHOMCH.'KNKOUS  MFD l A 


In  the  previous  section  a discussion  was  presented  of  two  poly- 
nomial models  for  transmittance  through  homogeneous  media.  Since 
the  atmosphere  is  an  Inhomogeneous  media  It  Is  necessary  to  somehow 
modify  the  previous  theory  to  make  the  results  applicable.  Only 
under  very  special  and  impracticable  assumptions  can  a model  he 
derived  for  the  inhomogeneous  case.  Several  procedures  are  discussed 
in  this  section  for  effecting  the  transition  from  inhomogeneous  to 
homogeneous  cond i t ions  . 


3.1  The  Transmittance  Function 

For  paths  through  the  Karth's  atmosphere  use  is  commonly  made 
by  the  "hydrostatic  approximation"  to  convert  the  integration  in 
(2.1)  to  an  integration  over  pressure.  Such  approximation  is 

du  - dr  (3.1) 

S 

where  M is  the  mixing  ratio  of  the  absorber  and  g is  the  gravity 
acceleration.  Beer's  law  for  an  inhomogeneous  path  from  P.  to  P„ 

l l. 

then  becomes 


Tv(P,T,u)  ■*  e 


k (p , t)m(p) dr 

V 


(3.2) 


which  must  also  be  averaged  over  a spectral  band  as  in  (2.3). 

In  view  of  essential  difficulties  in  the  exact  evaluation  of 
(3.2)  for  a path  along  which  P,,  is  significantly  different  from  Pj , 


■ - - 


use  must  bo  made  of  t ho  knowledge  gain  from  Its  evnlviat  ton  lor  homo- 


geneous paths. 


3.2  Method  of  Welureh  and  Neuendorf  I ei 

A given  atmosphere  may  he  described  by  a family  el  N curves- 
of-growth  representing  the  transmittance  along  a vertical  path.  Kach 
curve  is  for  constant  presure  P ^ and  temperature  T^,  where  i - 1,..., 
N,  as  depicted  qualitatively  in  Fig.  3.1.  Since  the  transmittance 


along  the  path  to  a given  pressure  level  I must  be  unique,  then 


T(pfT,,u;_i)  - riP^.T^.  Ul_j) 


t 3 . 3) 


where  n^  ^ Is  the  absorber  amount  contained  between  the  atmospheric 
top  and  level  l-l  at  conditions  of  level  1-1,  and  u‘  ^ is  the  same 
variable  but  at  conditions  of  level  i.  The  amount  uj  | Is  said  to 
he  "equivalent"  to  u‘  ^ in  the  sense  of  satisfying  ('.3).  According 
to  the  method  proposed  bv  Weinr«'b  and  Neuendorf  t or  ['J|  a solution  is 
first  found  for  u^'  ^ from  (.!•!!  with  both  a knowledge  of 
u ^ ),  and  an  available  band  model  for  homogeneous  paths  evaluated  at 
Pj  and  Tj.  It  follows  then  that  the  transmittance  to  the  next  pres- 


sure level  Is  obtained  as 


O'YW  - 1 < I , ’ T , + A,.p 


1 1. 


where  An ^ is  the  incremental  absorber  amount  contained  In  the  path 
length  AZ ^ between  the  pressure  levels.  This  incremental  gas  amount 
is  computed  with  a knowledge  of  the  absorber  density  p with  the 
equal  ion 


TRANSMITTANCE 


Au,  * p(T1,l’1)AZi 


(3.5) 


Although  la  principle'  (3.3)  is  exact  for  general  transmittance. 
Its  application  to  actual  cases  Involves  the  use  of  both  an  existing 
band  model  t(1’,T,u)  for  homogeneous  paths,  as  well  as  of  a solution 
for  at  each  level.  In  addition  to  intrudin'  t i ng  inaccuracies, 

these  steps  impose  additional  demands  on  the  computational  algorithm 
for  transmittance. 

3.3  Method  of  Curtis  and  Godson 

A classical  and  widely  used  alternative  to  the  above  method  is 
that  which  uses  the  Curt i s-Godson  [10,11]  relations  in  the  computation 
of  the  equivalent  gas  amount  u‘_j.  This  method  establishes  two  equa- 
tions tor  (3.3)  which  are  exact  only  In  the  limiting  conditions  of 
weak-and  strong-line  transmittance.  From  these  follow  equivalent 
relations  for  the  gas  amount  and  as  well  as  for  the  level  pressure. 
That  is 

T(ri*Ti*‘,i)  “ \o(Tl'UI-l) 

and 

Til’j.Vj.u^  « ^(''j'Tj.u  j_j)  (3.7) 

which  lead  directly  to  a solution  for  u^'  from  (3,b),  and  to  a 
solution  for  the  equivalent  level  pressure  r ‘ I t orn  ().7).  Such 
solutions  are  independent  of  the  type  ot  band  model.  They  follow 
directly  from  the  integral  form  of  Beer’s  law  (2.1)  and  the  l.orcntzinn 


half-width  (2.11)  as 


i-1  S, 


S (z)du(z) 


(3.8) 


1 Siui-l 


S(z)a(z)du(z) 


(3.9) 


whore  the  integration  is  carried  over  the  path  length,  and 


PA/To\l/2 


'o  h 


(3.10) 


The  transmittance  to  the  next  pressure  level  follows  according 


to  this  method  as 


T(Pi,Ti,Ui)  “ T^Pi,Ti’Ui-l  + ^ui"* 


(3.11) 


Briefly  stated, (3.11)  is  based  on  the  determination  of  an  equivalent  gas 


amount  following  a curve-of-growth  at  an  equivalent  pressure,  but  at 


the  same  temperature  as  that  provided  by  the  Weinreb-Neuendorf fer  method. 


A little  more  conceptual  development  may  be  used  [12]  to  show  that  the 


methods  are  identical  when  xfP^.T^.u^)  is  of  either  purely  weak  or  purely 


strong-line  characteristics. 


3.4  Method  of  McMillin  and  Fleming 


A recent  approach  to  the  problem  by  McMillin  and  Fleming  [13]  led 


to  the  separation  of  (3.4)  into  a product  of  two  functions,  one  of  which 


is  given  by  (3.3).  The  other  is  a four-parameter  polynomial  f expressed 


in  terms  of  temperature  differences  within  isothermal  layers  of  uniformly- 


mixed  gases.  That  Is 


(3.1.1) 


l 1 


T(pt*Ti*ui)  “ 1(p1,Ti,ui-i)f(ATi’ATrATr> 


i(pi-i,Ti-i,ui-i)f(ATi,,'TrATr) 


(3.13) 


where  AT^  is  the  level  temperature  difference  from  a reference  average 
temperature  , while  AT ^ and  AT ^ are  these  same  temperature  differences 
but  weighted  by  the  level  pressure  and  averaged  over  the  entire  atmo- 
spheric region  above  the  level.  Since  thev  follow  from  the  use  of  the 
equivalent  mass  concept,  they  may  as  well  be  called  "equivalent  tem- 
perature differences".  Specifically, 


ATj  - Tt  - T (3.14) 

jAT(P)dP 

AT'  - J (3.15) 

| dP 


PAT(P)dP 

AT"  - ' (3.16) 

| PdP 

wliere  the  Integration  is  carried  along  the  path  ending  at  P.. 

Tb.e  polynomial  f in  (3.13)  is  a function  representing  the  trans- 
mittance through  toe  layer,  and  is  given  for  a layer  as 


f - b,  + b,AT  + b.AT“  + b. AT  * + b -AT  (3.17) 

where  the  bj's  (J»l,...,5)  are  spectral  parameters.  In  particular,  b^ 
represent  the  ratio  of  the  transmittance  to  the  bottom  of  the  laver  to 


u 

the  value  at  the  top  for  the  reference  temperature.  The  original 
Intent  in  arriving  at  f was  to  allow  for  a convenient  procedure  for 
calculating  the  transmittance  through  a layer  at  any  temperature  T , 
with  a knowledge  of  the  transmittance  at  a reference  temperature  T^. 

J.5  Method  of  Transmittance  Product 

The  separation  of  (3.4)  into  the  product  form  (3.12)  suggests 
another  separation  in  the  form 

l(ri,Ti,Ui)  “ T(1>1»Ti»u[_1)AT(P.,Ti,Au  ) (3.18) 

“ l(Pi-l,Ti-l,ui-l)AT(ri,Ti,Aui)  (3.1*)) 

where  At  is  the  transmittance  through  Au^.  Even  though  the  approximations 
inherent  in  (3. Id)  are  of  the  same  magnitude  as  those  in  (3.12),  the 
function  At  is  not  restricted  to  temperature  dependence  alone, and  may 
be  used  tor  nonunlformly-mixed  gases  as  well. 

Further  attractiveness  of  (3. Id)  lies  in  the  fact  that  Al  may  bo 
represented  with  any  general  band  model,  which  shows  the  explicit  de- 
pendence on  pressure,  temperature  and  gas  amount . One  set  of  spectral 
parameters  describing  sui h model  could  be  developed  for  the  entire 
atmosphere  for  each  spectral  band  of  interest.  An  examination  of  (3. Id) 
shows,  additionally,  that  it  is  exact  for  the  monochromatic  case  and 
that  it  does  not  require  the  use  of  anv  equt\  donees  in  the  variables 


i nvolved . 


BAL'D  MODEL  DEVELOPMENTS 


1 S 


By  model  development  is  usually  meant  the  determination  of  the 
spectral  parameters  of  the  model  using  measured  transmittance  data. 

In  the  absence  of  a sufficient  amount  of  data  i rom  experimental  analyses 
use  is  mostly  made  of  line-bv-line  calculations  for  theoretical  paths 
resembling  the  actual  atmospheric  situations.  By  line-by-line  calcu- 
lations is  meant  the  evaluation  of  Beer's  law  in  the  form  of  either 
(2.1)  or  (2.2)  using  the  parameters  of  each  significant  line,  and 
assuming  some  type  of  broadening  function.  Tills  procedure  is  discussed 
In  detail  in  the  second  part  of  the  final  report  under  this  contract. 

In  this  section  a least-squares  technique  is  applied  to  the  determination 
of  the  spectral  parameters  for  the  two  polynomial  models  presented  in 
Section  2. 


4.1  Spectral  Parameters  tor  Weak-Line  Function 

The  weak-line  approximation  for  the  Maver-Goody  model,  as  well  as 
for  most  models,  was  given  by  (2.6).  In  terms  of  the  variable  in  (2.12) 


it  becomes 


w 


expi-k^u' 


(4.1) 


The  curve  fitting  technique  to  be  used  in  the  determination  oi  k consists 


of  mtnimiaing  the  difference  function  d as 


d(k> 


(4.2) 


where  the  summation  is  over  all  levels  for  which  data  exists  on  P 
Ti  and  Ti>  and  W£  is  the  weighting  function  [14] 


W = Ti  Ti)b  (4.3) 

For  convenience  in  the  minimization  process  the  change  of  constants 

B = 1/k2 

may  be  introduced  so  that  the  least-squares  derivative  of  (4.2)  becomes 


(4.4) 


4.2  Spectral  Parameter  for  Strong-Line  Function 

The  strong-line  approximation  of  greatest  versatility  is  the  one 
given  bv  King  and  stated  in  (2.9).  In  terms  of  the  variable  in  (2.13), 
it  becomes 


i 

s 


(4.5) 


Since  the  spectral  parameters  n,  and  C appear  within  the  function  G in 
(4.5),  their  determination  requires  a non-linear  curve-fitting  procedure 


The  quantity  to  be  minimized  Is 


where 


(4.6: 


Wi  “ 1 


i 


i , 
l 


.x 

s 


,1 


1 


' S ,i 


“ G(ntC) 


U iS  first  assumed  il'  <4-b>  th-'t  an  assumed  point  n\C*  is  very  clos 

to  the  point  which  yields  the  computed  minimum  of  d and,  then  d is 

* i 

expanded  in  a two  term  Taylor  series  about  the  assumed  point.  That 
is , if 


B » £n  C 

d (n,Bl  = Cn  a - (>\  * 

1 t s , i 


(4.7) 


then 


3d , 


dt(n.B)  •=  d.  (n\B  ')  + i 


(n-n  •)  + 


n * , B ‘ 


?B 


(B-B  ') 


n",B 


(4.8) 


After  substitution  of  (4.8)  into  (4.6), 


the  minimum  point  is  obtained 


through  the  derlvat iv  fS 
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(n,B)  « 0 


which  yjild  the  linear  sot 


#Zx* + ‘Z1!1” <>s-1  ■ i”a<} 


(A. 9) 


£vi  + A"V}i  |ln  o8fi  (n'.B' 

i i i ^ 


)~  In  >1  (4.10) 


In  this  set  of  equations,  the  following  notation  is  used 


x . = ~ — la  . 
v ^n  s,  i 


n',  B ' 


TT  In  a 
3n  s,i 


n , B 


An  *»  n - n * 


AB  = B - B 


In  matrix  form  (4.Q)  and  (4.10)  may  be  written  conveniently  as 


lp  I (A]  - [nl 


(4.11) 


with  solution 


[A]  - [P-il  [Dj 


(4.12) 


whore 
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4.3  Spectral  Parameters  for  Polynomial  in  Pressure,  Temperature  and 
Concentration 

In  a previous  section  a polynomial  model  was  derived  which  ex- 
pressed transmittance  as  a function  of  the  weak-and  strong-line  approx- 
imations. Although  the  model  as  given  by  (2.10)  allows  for  an  expan- 
sion in  products  of  functions  of  t , t , it  is  seldom  advisable  to 

w s 

include  more  than  the  first  three  terms.  This  premise  is  based  on 
the  fact  that  for  purely  weak-or  purely  strong-line  data,  the  models 


for  these  limits  normally  curve-fit  sufficiently  well  to  the  data. 

The  problem  lies  in  the  region  in  between  these  limiting  condition, 
since  the  available  data  r is  more  than  likely  an  inseparable  mixture 
of  weak-and  strong-line  transmittances . Any  additional  term  beyond 
the  third  will  tend  to  place  undue  weights  on  t lie  limiting  functions, 
decreasing  the  effect  of  the  "interpolating"  third  term. 

In  a three-term  expansion,  (2.10)  may  be  re-written  in  the  form 

2 2 

**  P x + py  + B XV  (4.13) 

w s w , s 


which  yields  the  least-squares  relation 


who  ro 


MM 


£V:  £v;V;  £v,V 


i i n i 


X •»  ■■—■ 

I 1m  , 
w,  I 


J l lm 


W * t ’.  (Ini  ,V 
f i.  t 


As  hoi oro . t ho  so  In i ion  lollows  ns 


I A | ■=  [11] 


(4.1 b) 


Spootr.il  r.it.imoto rs  lor  l'olynomi.il  in  Pressure,  Tempo  mt  tiro  ntul 
Concent  ration 


Tho  model  of  Smitli  expresses  t r.tnsmi  t tattoo  .is  .i  roivnomi.il  function 
of  fho  K.is  amount,  tho  pressure  and  tho  tompornturo . in  its  most 


usolnl  form  1 9 I tho  1’olvnominl  is  expanded  into  14  terms  ns 


» * 


111  (-lilt)  ” .1,  + it  X + « v -»  a./  + it_  XV  -f  .1  xz 
I 2 ' 4 5 f, 


2 2 32 

+ a x + a ux  z + a yz  + a,_x  + a.,xz  f a.  ,z 


7 8‘ 


+ .t,  .xv/  + a,  , vx‘ 
l J I a" 


10 


11 


(A.  16) 


For  ;t  least-squares  analysis  of  (.  4 . 1 6 > a minimization  Is  made  of 
the  difference 

‘Ha,.  j-*!,...,  14)  - ) ^ln(-lnr)  -g(a,,  |-1 14)^  (4.17) 

1=1 

where  g is  the  right-hand  side  of  (4.16).  Hoeause  ol  the  difficulties 
in  arriving  at  a weighting  1 unction  of  physical  significance,  no 
weighting  function  is  used  in  (4.17).  In  view  ol  the  sizes  of  the 
matrices  assoc iated  with  a solution  ot  (4.17),  they  ate  not  presented 
in  this  section.  They  are  included  In  the  Appendix  with  the  associ- 
ated computer  programs.  The  matt  ices  follow  the  derivatives 


*S~-  ‘Ha.,  .1-1 14)  - 0 

<*«  . 


(4.1 8) 


Although  the  matrix  solution  (4.12)  is  applicable  to  the  least- 
squares  set  of  equations  obtainable  from  (4.17).  a different  approach 
may  he  taken.  For  each  pressure  level  an  equation  may  he  written 
having  the  form 


a,  + n,  x,  + a,  y + . . . f a.,  y.x, 


(4.14) 


i 


tj  **  ln(-lm  ) 


x ■ l n u j 


y{  n 1 n (V^  /l'^ ) 


*1  * ln(W 


Define 


[Y] 


[H| 


1 x,  y1  ...  y^J 


1 X r 


y2  ...  y,x2 


lx  V 

N * N 


>’nxn 


[A] 


a. 


so  that  we  mav  write 


IY]  - [B]  [A | 


(4.20 


The  application  ol  (4.18)  to  (4.17)  yields  a set  of  14  equations 


of  the  form 


fcl  " "l  > 1 + *2  \ *i  + •• 


■ «i\\ + + 


I • 


which  using  matrices  m.iy  ho  written  as 

[S|  - [U?  [VI 


whore 


— ^ 

al 

[Si  - 

i 

, [V]  = 

a2 

ZbV. 

i 

l— 

a14 

L - 

[>’) 


The  sclntionof  (4.20)  Is 

IV]  - UJ]"1  [s] 


To  find  the  soldi  ion  lor  [V)  by  matrix  operation  on  the 


(4.21) 


previously 


defined  matrices,  note  that 


.'(> 

Am  uwrioN  10  i s»im  oo,:  uomooknkoi’s  paths 

In  order  to  evaluate  the  validity  .nut  usefulness  o t t ho  formula 
t ion  presented  in  tho  earl ior  soot  ions  ot  this  report,  sevetal  ot  t ho 
procedures  and  models  wots-  developed  lot  computet  solut  iotts.  In  this 
soot  Ion  a disoussion  is  presented  ot  tho  equal  ions,  procoduios  and 
oomputov  programs  assooiatod  with  tho  dovo  I opmottt  s of  thi'  homogeneous 
path  cast'  tt'r  t h«'  l llun  iV , spootral  hand. 

s . 1 Transm i 1 1 anoo  Pa t a 

l'ho  transmittance  data  tisod  in  tho  development  t'l  tho  hand  models 
tot  homogeneous  paths  consisted  o!  .'00  values  in  tho  ragnos  t rom  O.OOo 
to  0.1  attd  O.s  t o ‘>il  atm.  cm.  at  Standard  Temperature  anti  Pressure. 

Tho  data  wort'  available  1 t>i  all  six  ohatutt'ls  t'l  t Ito  vortioal  tempera 
tore  profile  retrieval  (VTPR)  experiment  in  the  llpm  TO,  hand.  fht' 
eenter  frequencies  tor  those  channels  .trt'  approximately  ('(>',  i'7‘,  i'94, 

708,  721,  and  747  cm  *.  Petal  Is  t't  the  1 ine-hy  -line  ealoulat  ional 
methods  ttsed  to  arrive  at  tht'se  ilata  are  given  In  part  two  ot  this 
final  report. 

1.2  Pcve lopment  of  Model  in  Fq.  4.11 

The  hand  models  chosen  for  curve-fitting  t o the  transmittance 
ilata  were  those  of  the  polynomial-  t vpo  as  proposed  hv  Pierlnlsst 
(Kq.  4.1.1)  and  Smith  (Fq.  4.1h).  In  particular,  4.11  had  hoen  proposed 
and  developed  some  time  prior  to  the  period  covered  hv  litis  contract 
[ 1 1 1 and  called  the  Five -Parameter  model,  The  Theory  ol  its  develop- 
ment is  repeated  here  tot  convenience  and,  at  the  same  time,  some  honol 
fictal  changes  to  the  original  version  are  introduced.  The  innovations 

,1 


( 


n 


are  incorporated  into  the  original  program  called  IRABSMD,  and  the 
resulting  program  renamed  ElkClM. 

In  the  earlier  work  the  exponential  limiting  function  for  hail 
been  substituted  into  the  comple  model,  in  such  a way  that  its 
parameter  K was  part  of  the  quadratic  parameter  11  . This  precluded 
the  curve-fitting  of  weak-and  strong-line  models  independently,  followed 
by  a curved-fitting  to  the  general  quadratic  model.  In  principle,  this 
latter  procedure  is  more  nearly  correct  or,  at  Least,  more  consistent. 
Initially,  only  the  strong-lino  model  was  curve-fitted  individually. 

The  following  is  a summary  of  the  model  version  developed  under  the 
present  effort . 

The  model  forms  developed  following  the  procedures  of  Section  4 

are : 


Weak-Line  (Determination  of  k) 


1 * exp 

w 


{-k  S u) 


(4.1) 


Strong-Line  (Determination  of  n and  C) 
t 


1 /n  ] 


i - l-G  { n,  [nr  (n)  P CS  (P/P  ) (T  /T)1/2 
a v.  L L tt  o o 


1/2 


f ( s . 5 ) 


Complete  Model  (Determination  of  B , R , and  11  ) 

w 8 ws 


f^-V 

In  I I 


- B 


-1 
In  t 


+ B 


In  T 


+ B 


ws 


In  T In  T 


l 4. 13) 


* ' • t s si 

Because  of  the  presence  of  the  incomplete  gamma  function  the  mo- 
del has  proven  to  be  somewhat  harder  to  develop  than  the  other  polyno- 
mial-type model  discussed  below.  The  problem  arises  in  reaching  con- 
vengence  for  the  strong-line  parameters  n and  C.  This  model  is 


HMKH 


non-linear  in  M and  C and  had  to  bo  linearised  for  its  present  appliea 
tion.  Thus,  the  determination  of  these  spectral  parameters  need  to  be 
determined  through  iterations  from  an  initial  guess.  This  is  a trouble 
some  and  time-consuming  process  which  is  required  every  time  that  the 
data  is  changed  and  new  parameters  need  to  he  developed.  For  this 
reason,  the  presentation  of  its  development  here  is  not  intended  lot- 
use  by  White  Sands  Missile  Range  in  their  V ITR  experiment.  For  that 
purpose  the  model  of  Smith  is  recommended,  together  with  some  technique 
for  inhomogeneous-to-homogeneous  path  conversion.  The  model  was, 
however,  developed  and  the  results  for  channel  1 is  included  in  Appendix 
A together  with  program  KIKl'lM.  The  results  have  not  been  optimised. 

5.3  Development  of  Model  in  Eq . -*.lb 

The  model  of  Smith,  discussed  in  detail  in  sections  3 and  4 ot 
this  report,  was  selected  for  development  using  the  homogeneous  path 
data  of  Sub-section  5.1.  This  model  is  easier  to  use  because  it  is 
simply  a linear  combination  of  logarithmic  functions  involving  the 
variables  u,  P and  T.  Calculations  using  the  developed  model  are 
equally  fast  and  accurate  since  the  process  consists  of  the  evalua- 
tion of  a polynomial  of  the  third  power. 

The  form  of  the  Smith's  mode!  used  is  given  by  4.1b  with  the  terms 
propsed  by  Vetnreb  and  Neuendorf f rr  l*1  1.  that  is 

2 2 

n(-ln  t)  - a,  + a„x  + a„v  + a.z  + a.xv  + a , x r + a,x  + acx  e 
id  3'  s b e b 

3 2 2 2 

a_vz  + a,„x  + n,.x2~  a,  . z " -♦  a,  ,xvr.  + a,  vx  fs.lb'' 

9 10  11  l-  13  14 

where  now 

x » 0.1  lit  u 
v - In  V 
2 - In  T 
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Equation  4.16  was  used  in  program  EIGGAM  2 for  the  determination 
of  the  spectral  parameters  a.  (i  « 1,  . . 14).  A copy  of  the 

program  is  included  in  Appendix  B together  with  a sample  of  the  out- 
put. The  results  for  the  six  VTPR  channel  and  for  all  principal 
absorbers  are  presented  in  the  next  Sub-section. 

5.4  Tabulation  of  Results 


The  transmittance  data  available  for  the  development  of  the  Smith's 
band  model  for  COj  and  H.,0  contained  tvco  ranges  in  the  gas  amounts, 
namely:  100  points  for  amounts  from  0.005  to  0.5  atm.  cm.  and  100 
points  for  amounts  from  0.5  to  50  atm.  This  break  down  was  selected 
in  order  to  cover  the  entire  portion  of  the  curve-of-growth  of  interest 
with  a high  degree  of  accuracy.  The  region  of  the  curve-of-growth  was 
decided  upon  on  the  basis  of  available  transmittance  data  for  a Standard 
vertical  profile  from  satellite  to  ground.  The  model  developed  for  the. 
latter  range  would  be  used  in  atmospheric  transmittance  calculations  to 
a larger  extent  that  the  former.  For  0^  the  range  from  0.005  to  0.5 
atm.  cm.  in  more  than  sufficient  in  order  to  be  able  to  calculate  the 
atmospheric  transmittance  for  any  reasonable  0^  profile. 

Tables  5.1  through  5.5  contain  the  14  Smith's  polynomial  para- 
meters for  all  ranges,  all  channels  and  all  gases.  For  calculations, 
these  should  be  substituted  in  Eq.  4.16  with  u in  atm.  cm.  t p in 
millibars  and  T in  degrees  Kelvin.  For  a given  homogenoeus  path  at 
these  conditions,  the  resultant  transmittance  will  be  the  product  of 
the  individually  calculated  transmittances  for  CO^,  ar>d  0^* 


r: 


M 


Once  determined,  the  spectral  parameters  of  Tables  5.1  through 
5.5  were  used  to  recalculate  the  original  transmittance  data  in  order 


to  have  some  measure  ot  their  correctness.  The  mean  of  the  absolute 
differences  between  the  calculated  and  the  original  transmi t tances 
were  evaluated  and  tabulated  in  Tables  5.6  through  5.11.  Included  in 
these  tables  are  also  the  peak  deviations  obtained  for  all  cases  con- 
sidered . 


fc.  . 


TABLE  5.1  SHITH  POLYNOMIAL  COEFFICIENTS  AT  STP  FOR  VTPK 

CHANNELS:  CC  IN  THE  RANGE  0.005  TO  0.5  ATM.  CM. 


TABLE  5.2  SMITH  POLYNOMIAL  COEFFICIENTS  AT  STP  FOR  VTPR 
CHANNELS:  CO,  IN  THE  RANGE  0.5  TO  50  ATM.  CM 
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Table  ">.12  Absolute  transmittance  dev  lit  I Ions  of  Smith  model  from 
original  data  using  Transmittance-Product  and  Curtis-Godson  method 
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TABLE  5.13.  SMITH  POLYNOMIAL  COEFFICIENT  FOR 


INHOMOGENEOUS  PATHS 


6.  APPLICATION  TO  CO.,: 

Previous  sections  Have  included  discussions  of  homogeneous-path 
models  for  calculating  the  mean  transmittance  over  a band.  For  the 
application  of  such  models  to  the  VTPR  experiment  it  is  necessary  to 
have  a means  for  conversion  to  the  inhomogeneous-path  cases.  Section 
3 covered  the  theory  ot  such  methods,  and  in  this  section  some  appli- 
cations are  discussed  and  evaluated. 

6.1  Transmittance  Data 

The  data  used  in  connection  with  the  development  of  models  for 
inhomogeneous  paths  consisted  of  100  transmittance  values  correspond- 
ing to  line-by-line  calculations  for  paths  from  the  satellite  height 
to  100  pressure  levels.  The  U.S.  American  Standard  Atmosphere  was 
used  in  the  calculations.  The  data  was  computed  by  NOAA  and  supplied 
by  the  Atmospheric  Science  Laboratory  at  White  Sands  Missile  Range. 
Only  the  CO,  absorber  data  for  the  six  V'lTR  Channels  were  made  avail- 
able for  model  developments.  The  gas  amounts  were  computed  using 
equations  which  were  suitable  to  the  particular  method  applied  in  the 
convers ion  t oehnique . 

6.2  Development  of  Cur t i s-Godson  Method 

Band  models  for  homogeneous  paths  mav  be  used  for  slant-paths  i! 
the  variables  u,  P,  T are  replaced  with  equivalent  variables.  One  ot 
those  equivalence  relations  that  are  available  is  t lie  one  provided  bv 
the  Curtls-Godson  relations  discussed  in  Section  1.1.  The  form  o. 
these  equations  that  were  used  for  t ho  absorber  amount  Is  obtained 


from  equations  0.2)  and  (3.8)  as 
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For  the  equivalent  pressure  (3.9)  and  (3.10),  that  i; 


i S,  u‘ 
t i-1 


S(z)  d(z)  du (z) 


1/2 


lead  to 


r T-'« 

i ,Siui/  n n n 


Au 


n 1 


41 


(3.8) 


(3.2) 


(6.1) 


(3.9) 


(3.10) 


(6.2) 


The  band  model  of  Smith  was  then  developed  using  the  slant  path  trans- 
mlttances  to  the  100  levels,  together  with  the  equivalent-homogeneous 
variables  u^,  P and  T..  The  identity  of  the  homogeneous-path  and 


inhomogeneous-path  transmit tances  was  established  in  (3.11). 


Program  G1GGAM  was  written  for  the  model  development  using  (b.l) 
and  (b.2).  This  program  is  essentially  the  earlier  presented  program 
EIGGAM  2,  but  mod i f led  to  account  for  the  calculation  of  the  equivalent 
quantities  required  bv  the  inhomogeneous-path  transmittance  data. 

Both  programs  appear  in  the  Appendix  to  this  report.  The  results  of 
the  calculations  are  presented  in  Table  5.12.  Appendix  D contains 
the  transmittance  data  used  for  the  inhomogeneous -path  model  develop- 
ment . 
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b.3  Development  of  Transmittance-Product  Method 

In  addition  to  the  method  of  Curtis  and  Godson,  three  other  methods 
were  discussed  in  Section  3,  namely:  the  methods  of  Keireb  and 

I 

Nuendorffer,  of  McMillin  and  Fleming  and  of  transmittance  product. 

For  use  in  connection  with  band  models  of  the  polynomial  type  the 
method  of  Woinroh  and  Neuendorf fer  is  difficult  to  apply.  This  is 
based  on  the  requirement  of  the  extraction  of  the  gas  amount  from 
the  avalytical  expression  for  an  existing  homogeneous-path  band 
model.  In  the  case  of  polynomials, this  involves  the  numerical 
solution  for  the  root  of  a tranc#dental  equation  at  each  pressure 
level  along  the  path.  In  addition  to  this  limitation,  the  band 
model  must  be  derived  from  homogeneous-path  transmittance  data  for 
condition  approximating  those  found  along  the  inhomogeneous  paths. 

That  is,  line-by-line  transmittances  for  the  inhomogeneous  paths  are 
never  used  in  the  development  of  the  band  model.  On  the  other  hand, 
the  method  of  McMillin  and  Fleming  is  restricted  to  uniformly-mixed 
absorbers  and  the  homogeneous-path  band  model  needs  to  be  developed 
for  each  pressure  level.  This  involves  an  enormous  number  of  spectral 
parameters  in  order  to  cover  the  entire  atmosphere  for  all  absorbers 
and  all  channels. 

The  method  of  transmittance  product  was  discussed  in  Sub -section 
3.5  and  consists  of  the  development  of  a band  model  for  the  incremen- 
tal transmittances  through  the  layers.  These  transmittances  are  given 
by  the  ratios  t(P{,  T{,  u.)  / rCP^,  T { t , vi  t ^ ) . as  given  by  3.19 
in  the  form 


.1 


For  At  the  model  of  Smith  was  used,  with  At  replacing  the  transmittance 
x.  Note  that  this  method  avoids  the  vise  of  equivalent  quantities,  and 
uses  the  inhomogeneous-path  transmi ttanees  in  its  development. 

Program  K1GCAM  was  modified  to  develop  (3.17)  and  compare  the  re- 
sults with  those  obtained  from  the  use  of  the  Cur t is-Codson  procedures. 
The  mean  and  peak  transmittance  deviations  are  presented  in  Table  5.12 
for  the  six  VTPR  channels.  The  model  parameters  are  tabulated  in 
Table  5.13. 


7.  DISCUSSION  AND  CONCLUSIONS 


The  principal  goals  of  the  portion  of  the  contractual  work  reported 
here  was  to:  a.  analyze  several  types  of  promising  band  models  for 
homogenous  paths  and  select  the  one  most  useful  for  the  application  to 
the  VTl’R  instrumentation,  b.  analyze  several  types  of  promising  methods 
for  applying  t lie  homogeneous-path  model  to  inhomogeneous  paths  and 
select  the  most  adaptable  to  part  a.  above. 

Band  models  of  the  polynomial  type  were  selected  since  it  is 
possible  to  improve  on  their  accuracy  by  the  addition  of  terms,  and 
they  can  be  developed  by  linear  least-squares  procedures.  The  poly- 
nomials proposed  by  Pierluissi  and  by  Smith  were  chosen  for  the 
analysis.  Both  of  them  were  subjected  to  the  same  input  conditions 
of  data,  computer  methods  and  the  same  feasibility  analysis.  Since 
the  polynomial  of  Pierluissi  includes  as  special  cases  a wide  variety 
of  models,  the  conclusions  derived  in  this  section  also  apply  to 
those  ramifications. 

The  polynomial  of  Pierluissi  from  a conceptual  point  of  view  is 
broader  in  applicability  and  of  greater  physical  validity  than  perhaps 
any  other  general  band  model  in  existence.  It  has  been  shown  In  the 
lierature  [6]  to  apply  to  high  resolution  data  with  high  accuracy.  a 

j 

With  the  use  of  King's  model  [31  for  strong-line  absorptance  it  is 
applicalbe  to  gases  having  absorption  lines  with  spectral  distribution 

i 

ranging  from  random  to  regular.  With  the  proper  selection  of  the 
spectral  hand  parameters  It  may  be  used  to  generate  most  of  the  classi- 
cal models  as  well  as  their  limiting  wrnk-and  strong-line  approxi- 
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mat  Ions . This  model  had  been  derived  earlier,  and  rederived  here  with 
little  change  for  the  purpose  of  documentation  as  well  as  for  incorpo- 
rating computational  modifications.  The  computer  program  to  be  used 
for  the  development  of  the  model  has  been  called  EIKCIM  and  it  appears 
in  the  Appendix.  The  developments  presented  in  this  report  are  limited 
and  are  incorporated  only  to  help  the  user  in  understanding  the  model. 
From  a practical  point  ot  view  the  model  is  restricted  by  the  inclusion 
of  the  strong-line  function  which  consists  of  an  incomplete  gamma 
function.  This  function  is  nonlinear  in  the  spectral  parameters  and, 
consequent 1\  , the  developmental  procedures  involve  the  linearization 
of  the  least-squares  equation.  It  is  necessary  to  iterate  in  these 
procedures  to  arrive  at  an  optimal  set  of  parameters  from  an  original 
guess.  In  general,  this  part  of  the  procedures  is  time  consuming  for 
the  user  and  wasteful  of  computer  time.  A second  limitation  worthy 
of  criticism  is  the  fact  that  the  model  variables  ( i.  e.  absorber 
amount,  presure  and  temperature)  appear  in  the  weak-and  strong-line 
function  rather  than  separated,  as  in  Smith's  model.  Tliis  forces  the 
user  to  go  through  three  developmental  procedures,  namely:  for  the 
weak-line,  strong-line  and  the  complete  model.  It  is  possible  to  do 
away  with  the  first,  but  not  with  the  second.  For  these  reasons  above 
it  was  decided  upon  the  recor.mendat  ion  and  usage  of  the  model  ot  Smith. 

Tlie  polynomial  model  of  Smith  was  programmed  as  program  EIGCAM 
(See  Appendix)  and  was  the  model  selected  for  the  application  to  the 
inhomogeneous  path  case.  The  model  was  developed  for  CO,,,  H^O  in  two 
ranges  of  gas  amounts,  and  for  0^  in  one  range.  The  average  transmit- 


deviations  obtained  for  the  six  V'l'PR  channels  over  all  the  ranges  for 


CO^,  Ho0  and  0^  are,  respectively,  0.00036,  0.00047  and  0.00017.  De- 
tails of  these  are  tabulated  in  Tables  5.6  through  5.11. 

With  respect  to  the  methods  of  conversion  for  homogeneous  paths 
to  actual  inhomogeneous  atmospheres,  the  results  indicated  that  the 
method  of  transmittance  product  should  be  used.  This  method  compared 
favorably  when  tested  together  with  the  method  of  Curtis  and  Godson. 

An  average  over  all  the  six  VTFR  channels  for  the  U.S.  Standard  1962 
atmospheric  profile  showed  a mean  transmittance  deviation  of  0.0004. 

This  compares  favorably  with  a corresponding  value  of  0.0011  obtained 
with  the  use  of  the  Curtis-Godson  equations.  Other  methods  were  analyzed 
in  Section  5,  but  not  implemented  due  both  to  the  scarcity  of  time  and 
to  the  fact  that  they  are  somewhat  impractical  or  inaccurate.  For 
instance,  the  method  of  equivalent  mass  is  only  a one-parameter  approxi- 
mation to  the  Curtis-Godson  method.  The  method  of  McMillin  and  Fleming 
involves  9,000  spectral  parameters  in  order  to  cover  the  six  channels, 
the  entire  atmospheric  and  the  three  principal  absorbers.  Finally,  the 
method  of  Weinreb  and  Neuendorffer  involves  the  development  of  the 
Smith  model  for  conditions  close  to  those  encountered  in  all  the  atmos- 
pheric layers,  as  well  as  the  iterative  solution  for  the  gas  amount  at 
each  pressure  level.  Although  the  results  of  the  effort  reported  here 
include  a band  model  for  uniform  paths,  it  was  developed  for  STP  condi- 
tions. In  spite  of  that,  the  iterative  solution  for  the  absorber  amount 
would  be  trotthlesome  to  obtain  because  the  Smith  polynomial  is  nonlinear 
in  that  variable. 
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APPENDIX  A 


1’ ROCK  AM  K1KCIM 

FIVE  PARAMETER  MODE1,  FOR  HORIZONTAL- 
PATH  TRANSMITTANCE  CALCULATIONS 

(Equation  4.13) 
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D G - 
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60  TO  ISO 

1 39 

0 N - 0 i 1 

6 0 ; o i s 3 

l 90 

IF  i J N . U E . 

60  TO  1 S 2 

1 s l 

o n - - o . : 

I S 2 

1> T,  ♦ 0 b 

N - N . i 

0 I - i)  H 

NT  -ON 


6 0 TO  V 

1 sO 

U = B * 0 » 

N = N + 0 N 

1 S S 

IF  ( N . 6 1 . 0 . 2 ) 
N-u  . 

60 

TO 

1 S 6 

1 *i  b 

it  i • b l • 6 • b ) 

N a 3 . ’ , 

60 

TO 

1 S 0 

1 o 

1C  0 

W H 1 T E ( 6 , 1 S 3 1 

) N 

• B 

1^31 

•c  0 f\  * 1 A 7 i 2 0 X | 3 H N 

- , F 

1 0 .s> 

0 0 1 b 2 K K — 1 , 

b 0 

1 -,7 

P 1 1 3 0 . J 

01=0.0 

0 0 l S o I a 1 , N I 

UlL  ABF{N,6*Al.06<US(I)  )/TL,PABLNFtPAOLBF  iAB) 

A - r*  A b L B F 
Pi:  PI  1 ♦ X • • 2 

IF  i Ab.Lt,.  0.0000001  1 (.0  fO  Isfl 
0 1 f - A U 0 o ( 1 . - T 0 ( 1 ) >/TL-AlO(s(  AO)  / T L 

0 1 - 0 1 ♦ X • D I F 
Iso  CONTINUE 

IF  l P l l . GT  .0 .0  ) 60  TO  l '19 
IU  b - 1 .0 

1 C » I C ♦ 1 

UilC.6T.ti)  60  TO  M 
6 0 T 0 1 S 7 

1^9  0 0 = 0 1 / p 1 1 

1 F < . . L E . 0 . b ) 6 0 TO  1 bO 

U 0 - 0 . i> 

uO  TO  1 5 1 

l«ju  : F l u ■ * 6 £ . ” 0 • 5 ) 60  To  1 b 1 
0u~-  « b 

1 b 1 b = b ♦ 0 B 

If  i -tiSiOfi  ) ,LT.  0.0005)  6 0 "o  SI 
1 5 2 CONTINUE 
b a - . 0 


; L A C - L I N E PARAMETER.  K 


!•;  i-o.o 
o i - u . o 

00  200  1 a 1 , N J 


I 


OIF--.  » / 1 A L Q G l T 0 l l ) ) * * 2 ) 

•Ml  .M  ; + v;  T * X * * *♦ 
01-D1'.jT»DIF»X«*2 
2 G 0 CONTI  NO  0 

c 

0 1 T ■ .M  i 

n = S ' T ( 0 £ T / 0 1 ) 

C 

c 

C iv.UA  VIC  P A H A M £ T £ K 6 S v,  , 0 S , B u S 

C 

C 

14  : 1 1 1 - 0 * 0 

P 1 2 = 0 • 0 
P 1 3 - 0 • G 
P 2 2 - 0 . 0 
P 2 j « 0 . 0 
P 3 3 -•  G • 0 
01  -0.0 
0 2 - 0 . G 
03—0.0 
M — 0 


\ITE16,14321N,B»K 

1432  .'OK'-  < 2 0 X 1 3 H N =tFlQ.6,SX,3N3  =>Fl0*S»5X»3HK  "*|j-  10.5) 

00  .K-i  I — 1 1 N I 

CAL-  A3FIN,B+AL0G(US( I ) )/TL,PABLnF,PABlBF,AB) 

1 r ( •L-.0.000CC01  •0F.A3.0E.  1 .OOOOUOO)  ^3  0 TO  16 

T = 1 .0 

. r i K • U .V  ( ! . G E . 2 5 • 0 ) GO  TO  16 
=-;./ALOG(EXP(-K*uW(I)  )) 
i — — . - / „ L C 0 t 1 . 0 - A 0 ) 

OIF-M./(ALOG(TO(I)  ) * » 2 ) 

; ’ 1 1 - p . 1 + v.  T * i X * * 4 ) 

P 1 2* P 12+  IV  T *<X**2MIY«*2) 

P l 3 - P 1 3 + ft  T • ( X * » 3 1 * Y 

P22-P22+WT  .(Y«*4) 

P 2 3 = P 2 3 + C T “ < Y 4 * 3 ) * X 


P3J  - P 3 J * T.  T 

• ( X**2) 

D 1 = D 1 - V.  T 

•0IF*X**2 

0 2 = 0 2 + A T 

•0  IF. Y«*2 

0 3 ■ 0 3 + '.V  T 

* 0 I F * X * Y 

- 0 N T . N U l. 

D-T-Ml*(P22*P33-P23*P23)-P12«‘iPl2«p33-P23*Pl3)+Pl3*(Pl2*P23-P22»P 

113) 

■3 - (Dl»lP22*P33-P23*P23)+02*<P23*Pl3-P12*p33)*D3*(Pl2*p23-P22 

; p 1 - ) ) / o £ T 

OS-  1D1*IP13»P23“P1.:*P33)+02»IP11»P33-P13*pl3)+D3«(pi3*pl2-Pli 

1 4 P 2 o > * / w L T 

d a <01*<P12*P23-\;2»Pl3)4D2MPl2»P13-Pll.p23).33*(Pll.P22-Pl 

. MP  1 2 > ) /OET 
c 0 TO  17  3 

171  DET=P1 1 •P22-P12«P12 

o.M(01«P22-D2*P12)/0£T 


tJi’  . i1 1 : - o ; • p 1 2 > / ol  ‘i 

ti  v « S 0 - 0 

Ni\-  . . - ; 

i . . ■ i_  V « *2  J 60  TO  1 F 3 

D l T “ ■ » 1 
3*'.  - . i/OcT 
6 b . . .0 

t>  ».  . U • w. 


T K . • M 1 T T A ,\  c £ TABLE 


A V l b - U . 0 
s V ..  . - U « 0 

•H  l - u < U 

PEA\S“G.Q 

w - ! - « 0 

PEAK  . i 0 


C = E X P ( 6 • T L ) 

00  . 7 ;*1,NU 

CALL  '\oF(N,6  + ALG6(US(  * ) )/”L,PAdLNF  ,PABLBF,A6) 

1 F ( ■> . • G.QGGGOOl  ) Ao-0*0l)UQ00» 

T S l I ) ; I • - A B 

T ..  I I ) - . . X P l - K • U iV  ( ; ) ) 

X = 1 . / { ,\  * u V.  ( I ) ) 

Y - - . . / A L 0 6 ( T b ( i ) ) 

LX  X« • 2 ♦Bb»  Y « • 2*fiL  S* \ * Y 


I F ( L 

X • u E • 0 • 

0 ) 60  T 

TC  ( I 

) ~E  XP  < - 

( 1 . / S Q R 

OTC  ( 

i ) =TC  ( 1 

1 -TO  ( I ) 

OTb  i 

, > = Tb  l I 

) -TO  ( I ) 

0 T V.  l 

i 1 -TV.’(  1 

) -TO  l 1 ) 

XML- 

N!  l 

AVLS-AVES*ABS(DTS<  I ) ) / X ' ; U 
AVlC-AVLC+A6S(07C(  • ) )/\..U 
A V i.  - A V t + A 3 S 1 D T IV  l 1 ) ) / A J 
I F < S(DTS(I)).GE.P;-AX3) 
i F 1 . 'j  ( 0 T C ( I ) ) • 6 E • P L A K C ) 

1 F ( . i S ( D T tV  ( I ) ) • 6 L • P E A K ;■  ) 


P£AnS  = A85(0TS<  I ) ) 
PEAXC  = ABS(i)TC(  1 ) ) 
P E A IUV  = A B * ( 0 T W ( 1 ) ) 


C 0 N T i N 0 E 

>J  R I \ t ..  0 i L , N , C , 6 ,V  , u S , 3 *•  S , X 

•'0.  ■ ■ i l H l , 3 H X , 1 2 H »V  A V L \ l ) ,\  o E R -,.7  ,/i31X,16hT£MPEKATUrE  = y C , / , 1 0 

. X , . ^STRONG  LlNt  PARAMETERS  , 1 2 X , 20 H QUADRAT  I C P A R A l T E R S » 2 X , 1 9 H >V  E A 
\ L . . 2 PARAMETER, /»1  IX,',  Fi  N . 3 * 2 x , 3 H C = , E l 0 • 3 , 2 X , 9 H >,  W a , ( 1 0 , 3 , 

2 X . R H S =,ElO»5,‘2X,bli.''b  =,E10«3,2X|3HK  -,ElG,b*//,2X,3n  LEvEl,  bX, 
• 1-  ' 0 0 E L VARIABLES, 18X.3YHCAR BON  0 1 OX  I DC  SATELLITE  T K A \ S M 1 T T A N C E 5 

/ X , 3 M N 0 • ,3X  , 9 H 1 1 L A N L1aE,2X,11hSTrON6  L 1 ME  , 7 X , 3 M 0 R I u I N A L , 2 X , 1 1 * S 
R i 1 L i N t , 2 X , H u \ V . u 1 2 x , * 0 H C A L C 0 L A T w 0 , 2 X * '•  1 1 0 \ , A T ' 0 , 2 x**  V H ft  E A 

• > L , 2X  , 9HDEV  1 AT  1 ON  1 t ) 

'o,  •.*)(!,  u,  t it,  , US  ( ; ) t * 0 i ^ ) , 1 v i ) , , 3 < 1 ) , i t ; ) , 

. 0 i ’ 1 , 7 2 ( I ) , D T >V  t I 1 , i • • ' , N U ) 

m j 1 \ 1 . t.  ,(i  • J 1 \ , F ° • H i X,F6*4i,bX,Fi#*i,iiX,Fo**,,*tX,F/ 


• 9 * ■ , r 6 1 *i  , ^ X , f , , 1 

a i •>2,PEANc,Pua\^,  u , A ^ l c t A ^ F ^ , A V i.  u 

0 a A 7 ( / / , 3 X , 1 0 h P 0 u Y X 0 M ; A L , l IX,  t l H S 7 N 0 N 6 LlNE»lSX,9m\£AK  L l N E , / , B X 


1 , l otiM  A a I MUM  EKlv'OR  , F8 . 4 , 3X,l3HMAXlMuM  E K K’  OK  , F d . H . s X , l 3 h (1  A X I n U M E R K 
••  0 i\  , , .S,/ibX*loHAVt.Ai>S«EH!K0R»F9«‘*»3Xil3HAVE*AaS»ESR0R,Fo«‘<,jX(ljH 
->  A V . A B S • £ K K 0 N i F a • H , / / ) 

CONTINUE  

STOP 

ENJ 


NO 


DIAGNOSTIC  • 
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SU  i,  ^OOTINt  ABF(N,X,PABlnFiPABLBF,AB) 

COi'l' ON/DATA/CAB  ( 19,  10) 

R A L N 

IF(X.GE»2.0#AND,N»Gt!»2,oO)  1 — 1 V 
I F ( \ . G i . 2 . 0 . A N l)  . N . L T . 2 , 0 0 » 1*18 

Ir(X.LT.2.0.AND.N.G£.2.00)  1=17 

IF(\.LT.2.G.ANC.N.LT.2.0U)  1=16 

1 ■ ( \ • L i *fc.*G«A\0*I'»*LT*l*2"5^  1 = 15 

IF  i X . L T . 1 •8.AN0.N.LT#  1 » 2 5 1 1 = 14 

I ( i . w.  T • . «O.ANO.N*GE.2,OG)  i = 1 3 
IFlX.LT.r.  0.AN0.N.LT.2.JU)  1 = 12 

1 . 1X.LT.t.U.ANG.N.LT»1.2‘j)  1 = 1 1 

l F ( \ . L T . 0 . 5 . A N U . N . L T • 1 . 2 5 ) 1 = 10 

lr  ( X * L T « 1 •0,AN0#N*GE#2,q0)  1 — 9 

IFIX.LT.0.0.AN0.N.LT.2.00)  1=8 

Ir*\.w1.0.O.AN0.N.LT«l«23)I=7 
1FIX.LT.-I..AN0.N.GE.2.0U)  1=6  , 

IF(\.LT.-l..AND.N.LT.2.00>  1=5 

IFiX.LT.-1..AN0.N.LT«1.25)  1=4 

IF(X.LT.-2..AN0.N.GE.2.0C)  1-3 

IF l X . L T . - 2 . . A N 0 . H • L T . 2 . 0 0 ) 1=2 

If  i\«Ll»*3«*AND<X*LT»l»  ‘bl  1 = 1 
TMl=2.»CA3(I,a)*N+2.«CAQ(I.4) 

T M 2 = C A B ( I , 2 ) ♦ C A B ( I , 5 ) * N + C A B ( I , 9 ) « N • N 
DIF=TM1**2-12.*CA3(  1 ,7)  * T M 2 
IFtDIF.LT. 0.0)  GO  To  1 
X .41  = ( - T M 1 ♦ S Q R T l D I.F  ) } / ( 6>  . * C A B l 1 ,7)  ) 

XM2M-TMl-SQRT(DIF)  j/(6.*CAB(  I ,7)  ) 

GO  TO  105 
1 X M l = - l 0 0 . 

XM2*  100.0 

1 0 3 X 2 = X » X 

ABLF  = CAB(I,l)-MCAB(I,2)t(CABU,5)+CA8(I,9i«N)*N)*X+ICAe(I,4)*CAB(I 
I ,7)  . X ) ♦ X 2 •»  ( C A B ( I , 3 1 ♦ C A B ( I ,8)*X2)*N«-(CAB<  I , o ) ♦ C A b < I , l 0 > • N > • N • N 

N 2 = 2 . * M 

PABlNF=CAB(  I , 3 ) ♦ ( C A B ( 1 , 5 ) ♦ C A B ( 1 ,8)*X)*X+(CAB(  I , 1 0 ) • i . 5 • n ♦ C A a t I ,6)* 
1 C A . ( 1 , 9 1 • X ) • N 2 
X 2 = 2 . • X 

P A 6 L B F = C A o ( I , 2 ) ♦ ( C A B ( I , 3 ) ♦ C A B ( 1 ,9>*N)*N*(CAS(  I , 7 > • 1 . 5 • X ♦ C A B < I 
IASI  1 , 8 ) * N ) • X 2 

IF  iASLF.oT.0.0.0R.X.GL.\M2)  GO  TO  2 0 8 
AB  =10.««ABLF 

If  ( as.lE.O.OOOOi  ) A8  = 0,o000i 
1 F l AB . GO . 0 . 99999 ) A8«0. 99999 

RETURN 

2 0 o AB  = 0 . 9 9999 

P A B L N F « 0 . 0 
PAbuSF ■ 0 « 0 
R ETUR N 
ENO 
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APPENDIX  B 
PROGRAM  EIGGAM  2 

SMITH  POLYNOMIAL  FOR  HORIZONTAL  PATHS 
(Equation  4.16.) 
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APPENDIX  C 
PROGRAM  EIGGAM 

DEVELOPMENT  OF  SMITH  POLYNOMIAL  FOR  INHOMOGENEOUS  PATHS 


c program  e icg am 

COMMON  Y ( 100  ) , £ ( li)u  ) » D ( 100  ) , NT  *WT(10(t) 

Oi'^SION  C(14)  *A(14f  15)  >IK(14)  »JC(14)  »UW(100)  rlJS(lOO)  »U(100)  » 

1 iAJ(io»i)  » S I ( 100)  »PI  ( 100)  #TI  (100)  * S ( 1 00  ) *P ( 100 ) * T ( 100 ) » TO 1 1 GO ) » 
X < 1 0 o ) 

L M 

Pu=l .01 325G+3 
T 0-273. 16 
NU=1 GO 


■ 


ooooo  r>  nor. 
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NF= o 

I”!-4*  • Rot\”4 
i>=VrtO  .»>? 

RmO=  i . *-kS2E—3 


input  Information 

KtAD(ft»  1 ) (PI  ( I ) » I=1#NU) 

1 FORMAT (10F7.2) 

WHlU(b,l)  (Pt(I)  ,I=1»NU> 
R£AO<h,?>  (T  f ( I ) , 1 = 1 » NU ) 
WHITE (n# 2)  (TI(I)»I  = 1»NU) 

2 FORMAT ( IUF6.3) 

DO  2S  J=1*NF 
WRITE (o» 222) 

222  FORMAT (1H1) 

R£rtO(b*3) (W#<SI<I>» 1=1,6)) 
WRITE  <r,,  3)  <w,  ( SI  < I)  » 1=1, 6)  ) 

3 FORMAT ( F10. 0»b£10.b) 

READ (5,4) ( T O ( I ) » 1 = 1 » NU) 
WRITER,**) (T0(I)*I=1»NU) 

*♦  FORMAT < SF1 0.7) 


ABSORBER  AMOUNTS 


a 


o 


7 


rt 


PUb=0.0 

PUTS=0.0 
DO  7 1=1, NU 


IF(I.GT.l)  GO  TO  ft 

TU)  = < <TI< I >+210.020)/ (T0*2.)  )**-0.5 
P(1 )=(PI ( I ) +3. 07E-4 ) / (P0+2. ) 

QUt I )=(M+(PI < I )-3.07l-4>*1.F+3)/(G*RH0) 


GO  TO  b 

T(I)=( (TI < I>+TI < I-l> )/(T0*2.) )+*-0.5 
P(1)=(PI ( I )+PI ( 1-1 ) >/<P0*2. ) 

DUl I ):|f>IPI ( I )-PI  ( 1-1 ) ) *1 .L+3) / ( 6*RH0 ) 
IFUI(I).GE.267.S)  S(I)=SI(1) 

IF ( T I ( I ) .GE.26P.5.AND.TI ( I ) .LT.287.5)  S(I )=SI (2) 
I F ( T I ( I ) .GE.237.ft.ANn.Tl ( I ) .LT.262.ft)  S< I )=SI (3) 
IF (TI ( I ) • GE .212.5. AMU . T I ( I ) .LT.237.ft)  S( I )=SI (4) 
1F(T  I { I ) .GE.187.ft.AND.TKl)  .LT.212.5)  S(I  >=SI  (5) 
IF ( T I ( I ) .LT • 1h7 . ft ) S ( I ) =SI ( 6 ) 

PHS=PUS+0U( I)*S(I) 


PU TS=PUTS+UU ( I ) *P ( I ) +T ( I ) *S ( I ) 
U.<(I)=PUS/S(I) 
US(I)=PUTS/(S(I)*P(I)*T(I>  ) 
CONTINUE 

t.R  ITE  ( b » 2 ) (P(I)  * 1 = 1 » Nil ) 


WRITE (6,2) <T( I ) »I=1»NU> 
WRITE (6,2) (S(I ) * 1 = 1, NU) 
V.KiT£(b,H)  (UW(  I ) , 1 = 1, NU) 
wRlTE<6,8)  (US ( I ) , 1=1 ,NU) 
WRITE(6»«)  (Ou(  I ) , 1 = 1, NU) 
Format ( iofio.5) 


nr,c  r>  r>  n n o o o r.  o 
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TRANSMITTANCE  DATA  SELECTION 


NT  = 0 

DO  0 I=l»NU 

IF(TOU)  .LT.0.0000001. OR. TO<  I)  .GT. 0.9999999)  60  TO  9 

NT-NT+1 

UF(NT)=UW(  I ) 

US < NT ) =U5 ( I ) 

□u < nt >=nu(  I > 

f 0 ( NT )=TO( I) 

WT (NT )=1 .0 
r (.\iT)=AL06(P(  I ) ) 

Z(NT)=AL06(T< I ) ) 

D ( i.T  ) r ALOG  ( ~AlOG  ( TO  ( NT ) ) ) 

9  CONTINUE 

POLYNOMIAL  coefficients  by  weak-line 


Du  °o  1=1# NT 
U<  i )=Ulv  < I ) 

99  x(i)=ALOG(U(I) )*0.1 
11=0 

10  Call  COEF(A.X) 

wRiTE(r>»  11)  ( (A(I.L)  »L=1»15)  #1=1.14) 

11  FGf<MAT(//#bX»l5Ffl.2> 

C‘*lL  LSIVEO( A, 14# IR.JC# 14# 1.0E-20#C#ERR> 

:i K 1 Tl  ( <>  # 12)  Err 

12  FORMAT ( bX  # lb ) 

.v  K 1 T F ( R » i3)(J#(C<I)»I  = l»l4)) 

14  FORMAT  (/  '/»25X#29HSMITH  POLYNOMIAL  COEFFIC  I ENTS  »//»3lX#  1 1HCHANNEL 
1 NO • * I4»/»2X# 2HC 1 # 6X  * PHC  2 » 6X » 2hC3  » 6X  * 2HC4 # 6X » 2HC 5 » 6X » 2HC6 . 6X » 2HC7 » 
2bX » 2UCR » oX • ?HC° » 6X  ♦ 3HC 1 0 # feX  # 3HC 1 1 » 6X » 3HC 12  # 6X » 3HC 1 3 * t>X  * 3HC 14  # / » 
Ol‘tL'J.4) 

TRANSMITTANCE  CALCULATIONS 
wkITF  (t»#  14)  J 

14  FOKMAT( 1H1»?5X»32HTRANSMITTANCE  FOR  CARHON  DIOXIDE »/#3bX» 11HCHANNE 
1L  nO, , 14  # / »?Xi 4HPATH#  3X  »6H AMOUNT » 3X » 8HPRESSURE » 3X  # 1 IhTEMPERATURE  » 
2UX.  1AuTRANSMITTANCE»/»3X#3HN0.  »3X»bHATM.CM*6X»2HM8#9X»lHK»4X»BH0F 
3 K>iNAL » 2X  » 10HC  ALCULATED  * 2X  * 9H0EVI  AT  I0N»  / ) 

C 

SCJ**i=0  • o 
DO  lb  1 = 1 » NT 

xP=C<l)+C(2>*X(I)+C(3)*Y(l)+C(4)*2<I)+C(5)*X(l)*Y(I)*Ct6)*X(I)*Z(] 

l)+C(7)*X(I)**^C(fl)^(I)*X(I)**2+C(9)*Y(I)*Z(I)+C(10)*X(ll**3+ 

2C(11)*X(I)*2(I>**2+C(12>*2<I)**2+C<13)*X(I)*Y(X)*2(I)*C<14)*Y(I>* 

3X(D**2 

IF  ( XP ,lE  .-2*> . 0 ) XP=-25.0 
IF ( XP. GT . +2b.  0 ) XP=+25.0 
XP=FXP ( VP ) 

IF(XP,lE.-2F».0)  XPS-25.0  I 
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IF(XP.GT.+2ft.O)  XP=+25.0 
TC=rlXP(-XP) 

u irc.Lf  .0. nun 01 ) rc=o.ooooo 

If (TC.GT. 0.99999)  fC=l. 00000 
nr=TC-TO( T ) 

SO'  I— SUM  + AGS  ( DT  ) 

.vRlTt  (<>•  lb)  I»U(I)»P(I)»T(I)»TO(I)»TC»DT 
lb  FGkMATOX#  I3»  1X*F9.4»  1X»F8.2»2X»FB.2»3X»F7.5»2X»F7.5»2X*F8.5) 

lo  CONTINUE 

AVLOT=bUM/FL0AT(NT) 

WRIT)  (o#l7)  AVEDT 

17  F OUV'AT  ( / / » rtX  # 19H  AVERAGE  DEVIATION  =#F8.5»//) 

11=11+1 

IF(II.GT.l)  GO  TO  1° 


POLYNOMIAL  COEFFICIENTS  BY  STRONG  LINE 


00  18  I = 1 » NT 

u< 1 >=uS< I > 

ltt  X(I)=ALOG(U(I) )*0.10 
GO  TO  10 
19  CONTINUE 


polynomial  coefficients  BY  transmittance  PRODUCT 


DU  21  I=1*NT 
U(D=DU(  I ) 

X( i )=ALOG(U< I ) >*0.10 
IF(I.GT.l)  GO  TO  20 
0( I >=ALOG(-ALOG(TO( I ) ) ) 

GO  TO  21 

20  DU)=ALOG(-ALOG(TO(I)/TO(I-l)  ) ) 

21  CONTINUE 

CALL  COEF ( A # X ) 

WRITE (6, 11) ( <A(I»L) *L=1»15) *I=1»14> 

CALL  LSIMEQ(A, 14» IR» JC»14»1.0E-20»C»ERR) 
WRITE (6. 12) ERR 

WRITE (A* 13) <J» (C( I) . 1=1*14) ) 


TRANSMITTANCE  CALCULATIONS 


L 


WRITE ( 6» 14 ) J 

SUM=0.0 

DO  24  1=1 » NT 

XP=C(1)+C(2)*X(I)+C(J)*Y(I)+C{4)*Z(I)+C«5)*X(I)*Y(I)+C(6)*X(I)*Z( 
l)+L(7)*X(I)**2+C(fl)*7(I)*X(I)**2+C(9)*Y(I)*Z(I)*C(10)*X(I)**3* 
2C(11)+X(I)*7(I)**2+C(12)*Z(I  >**2+C(13)*X(I)*Y(I)*Z(I)+CU4)*Y(  D* 

3X ( I ) **2 

IF(XP.LF.-2b.0)  XP=-25.0 
IFIXP.GF.+25.0)  XP=*25.0 
XP=tXP(XP) 

IF  < Xp.LE  »"2b. 0 ) XP=“25.0 


o n 
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IFUp.CE.+a^.O)  XP=+25.Q 
IFU.GT.l)  <30  TO  22 
TC=u'X?(-XP) 

GO  TO  23 

22  TC=TC*tXP(-XP) 

23  IF(TC.Lr.O.riOOOl)  TC  = 0. 00000 
IFtTC.GE. 0.99999)  TC=1. 00000 
3T=TC-T0< I ) 

SU.-.=SUr,+ABS(DT> 

.-/KITE  (6»  15)  I,U(I)  »P(I)  »T(I)  »TO(I)  »TC»OT 

2*.  CONTINUE 

AVEOT=SUM/FLOAT  < NT ) 

.•JR1TE(6»17)AV£DT 
2b  CONTINUE 
STOP 
END 


compilation:  no  diagnostics. 


SUBROUTINE  COEF(A#X) 

CON-’ON  Y ( 100)  ,2(100)  .0(100)  »NT»WT<  100) 
DIMENSION  A(14»15) »X(100) 

DO  1 1=1 .14 
DO  1 J=1 » IS 
1 A(ItJ)=0.0 

DO  2 1=1. NT 
WT(I)=1.0 

AU»1)=A(1.1  )+1.0*WT(I) 
AU»2)=A(1»2)+X<I)*WT(I> 
All»3)=A(l»3)+Y(I)*WT(I) 
A(1»4)=A(1»4)+Z(I)*WT(I) 

A(1#5)=A (l.S)+X( I > *Y( I ) *WT ( I ) 
A(l»6)=A(l»h)+X(I)*Z(I)*WT(I) 
A(1»7)=M1*7)+X(I)*X(I  ) *WT  ( I ) 
A(1.G»=A(1.H)+X(I)*X(I)*Z<I)*WT(I) 
A(1.‘>)=A(  1.BJ+YI  I ) *Z(  I ) *WT  ( I > 
A(1»10)=A(1.10)+X( I)»*3*WTII ) 

A ( 1 » 1 1 ) =A ( 1 » 1 1 ) +X  1 1 ) *Z ( I ) *Z ( I ) * WT ( I ) 

A ( 1 * 12 ) =A ( 1 » 12 ) +Z ( I ) *Z ( I ) *WT ( I ) 

A ( 1 » 1 3 ) =A { 1 » 13 ) +X ( I ) *Y ( I ) *Z ( I ) *WT ( I ) 
A(1»14)=A(1.14)+X(I)*X(I)*Y(I)*WT(I) 
A(1»15)=A(1»15)+D(1) *WT < I ) 

A ( 2 t « ) = A ( 2 . B ) +X  ( I ) **  3*Z  ( I ) * V T ( I ) 

A ( 2 » 10)=A(2* lQ)+X( I ) **4*WT ( I ) 
A(2»ll)=A(2»ll)+X(I)*X(I)*Z<I)*Z(I)*WT(I) 
A(2,i3)=A(2.13)+X(I)*X(I)*Y(I)4Z(I>*WT(I> 
A(2.14)=A(2. 14)+X( I )**3*Y(I)*WT( I) 
A(2.1S)=A(2.15)+X<I)*D(I)*WT(I) 

A < 3 . 3 ) =A { 3 » 3 ) +Y ( I ) *Y ( I ) *WT ( I ) 
A(3.5>=A(3»5)+X( I >*Y(I)*Y( I)*WT< I) 
A(3.S):M3.‘»+YII)*YII)*ZII)*WT(I) 


I | Mta 


A (3. 11 )rA(3» 11 l+Xli > *Y< I >*Z< I ) *7  ( I)*WT( I) 

A(3. 12>=A(3» 12>+Y< I > *Z< I >*Z< I ) +WT( I> 
A(3.13)=A(3»13)+X(I)*Y(I)*Y(I) *Z< I )*WT( I ) 

A ( 3 » 1 4 ) ( 3 • 14 ) +X ( I ) *X ( I ) *Y  < I ) +Y ( I ) ♦WT ( I ) 

A ( 3 » 15 ) = A ( 3 * 1 F» ) + Y ( I ) *0  ( I ) *WT  ( I ) 

A (4, li )=A(4»11 )+X( I ) *Z( I ) * * 3 * W T ( I ) 
A(4,12)=A<4#12>+Z<1>  * + 3 * W T ( I ) 
A(4»13)=A(4.13)+X(I)*Y(I)*Z(1)*Z(I)+WT(I) 

A ( 4 # 14 ) r A ( 4 » 14 ) +X  < I ) *X ( I ) * Y ( I ) *Z ( I ) ♦WT ( 1 ) 

A(4, 15)=A(4»  15)+Z  ( I ) *0(  I ) *.VT  ( I ) 

A(^»<t)=A(5»H)+X(  I ) ♦»3*Y  ( I ) *Z(  I ) *WT(  I ) 

AO>»10)=A<5»  10>+X<  I > **4*Y  ( 1 > *WT  ( I ) 
A<b.li)=A<b.ll)+X(I)*X(I)*Y<I)*7(I)*Z<I)**T<I) 

A (a* l3)=A(b. 13) +X ( I ) *X< I ) *Y(I ) *Y ( I ) *Z  ( I ) *WT ( I ) 

A(b»i4)=A(b.l4)+X(I)**3*Y<I)*Y(I)*WT(I) 

A(b,lb)=A(b.l5>+X<I>*Y(I)*D<I>*WT(I> 

A(o,H)=A(b»rt)+X(  I ) **3*Z  < I ) *7(  I ) *WT<  I ) 

A(o»m)=A  <6»«)  +X(  I ) *Y(  I ) *Z<  I >*Z  ( I ) *WT  ( I ) 
Alo(10)=A(o.ln)+X(I)**4+Z(I)*WT(I) 

A(o» 1 1 ) =A (o. 11 ) +X ( I ) *X ( I ) *Z ( I ) **3*rtT ( I ) 

A(u,l3)=A(o.l3)+X(I)*X(I)*Y(I)*Z(I)*Z(I)*WT(I) 

A(o»l4)=A(o#l4)+X(I)**3*Y(I)*Z(I)*WT(I) 

A(o»lb)=A(b.lb)+X«  I)*Z(I)*D(I)*tJT(I> 

A( 7, 10  >=A(7* 10 )+X ( I ) **5*WT ( I ) 
A<7,il)=A<7#ll>+X(I)**3*Z(I>*Z(I>*WT(I) 

A <7# 14)=A(7» 14 )+X (I ) **4*Y( I ) *WT ( I > 
A(7,lb)=A(7#lS)+X(I)*X<I>*D(I)*WT(I) 

A(6»8)=A (8 #8) +X ( I ) **4  + Z < I ) *7 ( I ) *WT ( I ) 
A(b,io)=A(«.lO)+X(I) **5*Z ( I )*WT( I ) 
A(b,U)=A<d.ll)+X(I)**3*Z(l)**3*WT(I) 
A(o,13)=A(6.13)+X(I)**3*Y(I)*Z(I)*Z<I)*WT(I) 

A(oi 14)=A(0» l4)+X ( 1 ) **4*Z ( I ) * Y ( I ) *WT ( I ) 
A(6,ib)=A(tt.lb)+X(I)*X(I)*Z(I)*D(I)*WT(I) 

A(9,9)=A (9,u)+Y( I )*Y< I )*Z( I)*Z<  I ) *WT  ( I ) 

A <9. 11 )=A<9. 11 )*X ( I ) *Y ( I )*Z( I ) **3*WT ( I ) 
A(9,l?)rA(9»l2>+Y( I )*7< I >**3*WT<  I) 
A(V#13)=A(y*l3)+X(I)+Y(I)*Y(I)*Z(I)*Z(I)*WT(I) 
A(9fl4)=A(9.l4)+X(I)*X(I)*Y(I)*Y(I)*Z(I)*WT(I) 

A (9,  lo)=A(9.  l‘j)+Y(  I )*Z(  I )*D(  I >*WT<  I ) 
A(10»10)=A(10,10)+X(I)**ft*tVT(I) 

A ( 1 0 » 1 1 )=A(J0,11)+X(I)**4*Z(I)*Z(I)*WT(I) 
A(10.14)=A(10.14)+X(I) **5*  Y ( I ) *WT ( I ) 
AUU*lb)=A(10,15)+X(I>*  + 3*L><I>*WT<I) 
A(llrll)=A(ll,H)+X(I)*X(I)*Z(I)**4*WT(I) 
A(11»12>=A(11,12)+X( I)*Z( I )**4*WT(I) 
A(11,13)=A(U»13)+X(I)*X(I)*Y(  I ) *Z  ( I ) **3*WT  < I ) 

A ( l 1 ,14)=A(11,14)+X< I)**3*Y( I ) *7 ( I ) *Z ( I ) *WT (I) 
A(llrlb)=A(ll,lb)+X(I)*Z(I)*Z(I)*D(I)*WT(I) 
A(12»12)=A(12*1?)+Z(I)**4*WT ( I ) 

A(12.13)=A(12, 13)+X(I)*Y( I )*Z( I)**3*WT(I) 
A(12,14)=A(12»14)+X(I)*X(I)*Y(I)*Z(I)*Z(I?  * W T C I > 

A(  12. 15)  =A  ( 12, 15)+ZC  I ) *Z(  I )*Il(  I ) *WT<  I > 

A(13.13)=A(13, 13)+X( I )»X( I )*Y( I )*Y( I)*Z( I)*Z(I)*WT(I ) 
Aa3.l4)=A(13»14)+X(I)**3*Y(I)*Y(I)*Z<n*WT<n 
A(13»lb)=A(13»15)+X(I)*Y(I)*Z(I)*D(I)*wT(I) 
A(14,14)=A(14,14)+X(I)**4*Y( I )*Y(I)*WT(I) 
A(14,15)=A(14,15)+X(I)*X(I)*Y(I)*D(I)*WT(I) 
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continue 


A (2 

-2>=A(1 

* 7 ) 

A (2 

-3)=A(1 

.5) 

A<2 

r 4 ) = A ( 1 

* 6 ) 

>\i£ 

bJ»A(l 

* 14) 

A ( 2 

6 ) = A < 1 

»«) 

A (2 

7)=A(1 

* 1 0 ) 

A (2 

9)=A< 1 

*13) 

A ( 2 

1 2 ) =A ( 1 » 1 1 ) 

A ( 3 

4 ) =A ( 1 

*9) 

A ( o 

6)=A(1 

*13) 

A ( 3 

7)=A(1 

*14) 

A ( 3 

S)=A (2 

*13) 

A (3 

10)=A(2»14) 

A ( 4 

4)=A(1 

.12) 

A ( 4 

b)=A(l 

*13) 

A (4 

6)=A(1 

*11) 

A ( 4 

7)=A(1 

*8) 

A ( 4 

<8)=A (2 

*11) 

A ( 4 

9 ) =A ( 3 

*12) 

A ( 4 

*10>=A(2*»> 

A (3 

r5)=A(3 

* 14) 

A ( b 

»6)=A (2 

* 13) 

A(b 

7)=A(2 

* 14) 

A ( b 

9)=A(3 

» 13) 

A ( b 

12)=A(4»13) 

A ( o 

6)=A(2 

#11) 

A ( o 

7)=A(2 

# 8 ) 

A ( 6 

1 2 ) =A ( 4 # 1 1 ) 

A ( 7 

7)=A(2 

*10) 

A ( 7 

8)=A(o 

*10) 

A ( 7 

9)=A(2 

*13) 

A ( 7 

12 ) =A ( 2 » 1 1 ) 

A ( 7 

13)=A<5#8) 

A(3 

9)=A (6 

*13) 

A(6 

12)=A( 

6*11) 

A ( 9 

10)=A(o*14) 

A(10»12)=A(7»11> 

A(10»13)=A(fl#14) 

DO  3 1=1 » 14 

JO  3 J=1 * 14 

A ( J » I ) =A ( I # J) 

RETURN 

END 


.ation: 


NO  DIAGNOSTICS. 


SAMPLE  OUTPUT 


77 


f 


: I 
*• 


TRANSMITTANCE  FOR  CARBON  DIOXIDE 
CHANNEL  NO.  1 


PATH 

AMOUNT 

PRES5URE 

temperature 

transmittance 

NO. 

ATM.  CM 

MB 

K 

ORIGINAL 

calculated  deviation 

1 

.0024 

.00 

1.16 

.99125 

.99126 

.00000 

• 2* 

<•  • UUbl 

.00 

1.20 

.98680 

,98880 

.00000 

3 

.0101 

.00 

1.1b 

.98553 

.96555 

.00003 

4 

.0206 

.00 

1.11 

.98093 

.96061 

-.00012 

5 

,031b 

.00 

1 .08 

.97542 

.97532 

-.00010 

6 

. 0492 

.00 

1.05 

.96916 

.96946 

.00028 

7 

,0o94 

.00 

1.03 

.96178 

.9ol97 

.00019 

a 

.0994 

.00 

1.02 

.95276 

.95274 

-.00001 

9 

.1348 

.00 

1.01 

.94166 

.94127 

-.00041 

10 

.1777 

.00 

1 .01 

.92645 

.92815 

-.00030 

11 

• .2262 

.00 

1.00 

.91306 

.91328 

.00022 

12 

.2914 

.00 

1.01 

.89552 

.89560 

.00008 

13 

, 3b46 

.00 

1.01 

.87594 

.87555 

-.00039 

14 

.4390 

.00 

1.02 

.85461 

.85425 

-.00036 

lb 

,5401 

.00 

1.03 

.83190 

.83160 

-.00029 

16 

«o563 

.00 

1.04 

.808 26 

.80830 

.00002 

17 

.7876 

.00 

1.04 

.78427 

.70486 

.00059 

IB 

.9367 

.00 

1.05 

.76025 

.70150 

.00125 

19 

1,1059 

.00 

1.06 

.73650 

.73618 

.00168 

20 

1.2954 

.00 

1.07 

.71314 

.71501 

.00188 

21 

1.4601 

,01 

1.07 

.69007 

.68941 

-.00067 

22 

1.6950 

.01 

1.06 

.66711 

.66740 

.00028 

23 

1.9552 

.01 

1.09 

.64400 

.64579 

.00179 

24 

2.2406 

.01 

1.09 

.62045 

.62321 

.00276 

2b 

2.5539 

.01 

1.09 

• 59b28 

.59039 

.00211 

26 

2.6949 

.01 

1.10 

.57136 

.57168 

.00032 

27 

3.2687 

.01 

1.10 

.54571 

.54403 

-.00168 

20 

3 • 6754 

.01 

1.10 

.51939 

.51596 

-.00342 

20 

4,1149 

.02 

1.10 

.49248 

.48796 

-.00452 

30 

4.5924 

.02 

1.10 

.46512 

.40003 

-.00509 

31 

5.1077 

.02 

1.11 

.43745 

.43208 

-.00537 

32 

5 , 6b09 

.02 

1.11 

.40961 

.40438 

-.00522 

33 

6,2570 

.02 

1.11 

.36177 

.37710 

-.00466 

34 

6.8936 

.03 

1.11 

.35409 

.35030 

-.00380 

3b 

7,5781 

.03 

1.11 

.32076 

.32411 

-.00265 

36 

6.305b 

.03 

1.11 

.29992 

.29851 

-.00141 

37 

9.0836 

.03 

1.12 

.27374 

.27378 

.00004 

3 a 

9.9122 

.04 

1.12 

.24636 

, 24999 

.00163 

30 

10.7937 

.04 

1.12 

.22393 

.22736 

.00343 

40 

11.7264 

.04 

1.12 

.20057 

.20580 

.00523 

41 

12.7166 

.05 

1.12 

.17638 

.16492 

.00654 

42 

13,7b94 

.05 

1 . 12 

.15747 

.10431 

.00665 

43 

14,8783 

. 06 

1 . 12 

.13788 

.14407 

.00619-a  . 

44 

lb.0504 

• 06 

1.12 

. 1 1 966 

.12463 

.00496 

4b 

17.2856 

.07 

l .12 

.10291 

. 10640 

.00349 

4b 

16,5890 

.07 

1.12 

.08760 

.08968 

.00208 

47 

19,9561 

.06 

1.12 

.07374 

.07463 

.00069 

40 

21,4005 

.08 

1.12 

.06135 

.Obi  34 

-.00001 

49 

22.913b 

.09 

1.12 

.05039 

.04978 

-.00061 

50 

24,5024 

.09 

1.12 

.04081 

.03987 

-.00093 

bl 

2b.lb7l 

.10 

1.12 

.03256 

.03150 

-.00106 

52 

27.912b 

.11 

1.12 

.02558 

.02451 

-.00107 

!■ 

K 

i 1 


I 


k 


53 

2y.73«9 

.11 

1.12 

.01975 

.01877 

-.00098 

54 

31.6465 

.12 

1.12 

.01495 

.01411 

-.00084 

55 

33.6441 

.13 

1.12 

.OHIO 

.01040 

-.00069 

5b 

35.7261 

.14 

1.12 

.00606 

.00751 

-.00055 

57 

37.9030 

.14 

1.12 

.00572 

.00531 

-.00040 

56 

40,1668 

.15 

1.12 

.00395 

.00370 

-.00025 

59 

42.5306 

.16 

1.12 

.00265 

.00255 

-.00010 

60 

44.9910 

.17 

1.12 

.00173 

.00173 

.00000 

61 

47.5498 

.18 

1.12 

.00109 

.00113 

.00004 

o2 

50.2123 

.19 

1.12 

.00067 

.00070 

.00003 

63 

52.9762 

.20 

1.12 

.00039 

.00040 

.00001 

64 

55,6503 

.21 

1.12 

.00021 

.00020 

-.00000 

65 

5b . 633b 

.22 

1.12 

.00009 

.00009 

-.00001 

66 

61.9279 

.24 

1.12 

.00002 

.00003 

.00001 

AVERAGE  DEVIATION  = 

.00156 

, 
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